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Abstract 


Let y = (p, q) be a canonical set of variables undergoing a 
near-identity transformation that involves a small parameter 


V 1 


z = 


.(k) 




(y) 


Then the transformation will be canonical provided v, 
has the form 


- (k) 


.(k) 


s 


C \ (k> 


f (k) (S) 


where v is a gradient operator in "conjugate phase space" 
y = (q, -p) , is arbitrary and may be regarded as the k-th 

order'of a generating function and is a given vector depen- 

ding on lower orders of £ and on their derivatives. Various choi- 
ces of f_ exist (all differing by some conjugate gradient) and 
some of these may be related to the conventional form or to Lie's 
form of canonical transformations. Properties and uses of this 
representation will be discussed. 


Equations are numbered according to the slide on which they 
appear in the presentation. 
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I would like to present here a somewhat new approach to the old 
subject of canonical transformations. 

I am sure everyone here is quite familiar with such transformations. 

If we start with a system of canonical variables 

y = (p, q) 

having N degrees of freedom (that is, y has 2N components) and 
seek a transformation to a new system 

2 = (£. 2) U-l) 

then the usual way of expressing this is by means of a generating function 
CT depending on N "new" and N "old" variables. There exist many forms 
for this, but we shall choose the form 

O’ = 0" (P, q) 

This gives the equations of transformation as 

= ^CT/'DP i 

( 1 - 2 ) 

p ± = o' / ^ q^ 

Now among all the various canonical transf omations , a particularly 
important type consists of the near-identity transformations, depending 
on a small parameter £ « 1 in a way which reduces them to identity 
transformations y=z. as £ -*• 0 . A generating function for such a 
transformation may have the form 


cr(l, q) = 


£ £ k cr (k! (P, q) 

(1-5) 

from which 

*i “ 

Y\k 

ii + /_, £ 

M< k) /DP i 

(1-4) 

Pi = 

\ ♦ tv 

^<T (k) /:>q i 

(1-5) 


This form is the basis of the familiar Poincare-Von Zeipel method: if 
we transform without time dependence, the new Hamiltonian H*(P, £) equals 
the old one H(p, q) , and by substituting the above formulas we get 



H 
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(p t + 2>W (k) /<V' qj 

= h* ^p ± » + 2,£ k ( d-<) 

This relation involves only the 2N variables (P, q) and by expan- /, \ 
ding, demanding that H* does not contain the angle variable and that (T' 
has no secular parts, solutions to problems in perturbed periodic motion 
may be worked out order by order. 

Unfortunately, this transformation from the original variables to more 
desirable ones is express** 4 , by functions of both old and new variables, and 
one still must untangle this relationship order by order. 

It would be much better if we had a direct transformation, of the form 


z = y + 2^ ^ k hy) 

(2-1) 

or, for many uses 


? = - * 2L £k T (k)( ^ 

(2-2) 


in which all new variables are explicitely given in terms of old ones, or 
vice versa. 

One way of getting such a transformation is by Lie*s method . About 
80 years ago Lie showed that if is a Poisson-bracket operator 

L^f) = [f, W] (2-3) 

then the transformation 

z. = exp ( e Ly) * y (2-4) 

(asterisk means operation) is a near-identity canonical one. The expo- 
nential is here defined by its Taylor expansion and may be reduced to a 
series like the one in equation (2-l) (or else, interchanging z and y 
gives a series like 2-2 ). If furthermore W can be expanded in powers" 

of £ 

w = w (k) (y) (2-5) 

then each order of the expansion will introduce the corresponding order of 
¥ , and will in addition contain assorted expressions involving lower orders. 



- 4 - 


This may then be sujslitvusd int.- .tie ,'ui notre-Vo... <Jeipel method, that is 
H(y) = H*(y + ? ^ k> ) 

and again separated into a set of expressions, one for each order. Just 
as ordinarily each order contains an undetermined function <T t * ,) which one 
adjusts to eliminate the angle variable and secular terms, so here one 
does it with V M . This is essentially the method advocated by Horl in 
1966 and later developed by Deprit. 

Sometimes, of course, one prefers to deal directly with the form 
given in equations (2-l) and (2-2) 

2 = y + (2-1) 

_y = z + (2-2) 

For instance, the Bogoliubov-Krylov method (or its mirror image used 
by Kruskal) leads to such a transformation; if y is a canonical set, 
one may then ask under what conditions is the transform., tion also 
canonical. 

I want to devote the rest of tha 

talk to this question : as you will see, it leads 

to a representation which includes both Lie transforms and the convent- 
ional generating function. We will deal with the form (2~l), but every- 
thing can be equally well done with the inverse series (2-2). 

Note first of all that eq. (2-l) deals with the vectors y and z_ 
without splitting them into canonical momenta and coordinates . ~We will 
try to preserve this unity throughout the calculation, and it can be done 
— but for a price: one must introduce "conjugate” vectors. 

Let us define them: if 

y a (p» q) 

then "y-conjugate” is 

y = (q. ~p) (3-i) 

That is, y is formed by rearranging the components of y so that 
every p^ ii replaced by its conjugate and every q^ by nanus its 
conjugate. 

With this notation, Hamilton's equations simply become 
dy^/dt * - DH/jy. 


(3-2) 
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or, in vector form 


dy/dt = - VH 


(3-5) 


Poisson brackets involve only one summation 



\'/ 't) a b ^ a \ 


V ~Db 

Z-^^yi 


Va-Vb (3-4) 


In particular 


and 


[»• 

[v y s] 


Q a/Dy^ 




ij 


(3-5) 

(3-6) 


which is the condition for a set y to be canonical. 

Now back to the main question* If y is canonical and z differs from 
it only by small terms, the condition for z, to be canonical is 


Lv z j ] - tv y j] - 0 


(4-1) 


Substituting 


|V £ k d k) - - IW 0 <«> 


. m v (m) 


This separates into a series of relations, one for each order. The zeroth 
order cancels, of course. The 0(& ) part gives 






or, by equation (3-5) 





w 


0 (4-3) 


0 (4-4) 


space. 


This looks very much like the curl dyadic of in y 

and in fact we shall designate it as such and write 

( t X = o 


(4-5) 


One can easily guess the solution: if the curl of a vector vanishes, it is 
the gradient of some scalar, so 


5 


(1) 


V X 


( 1 ) 


(4-6) 


with arbitrary. 


In general one finds, for the k-th order 

k-1 


(V x 


1J 


y (5-D 


m=l 


Since all terns on the right are of orders lower than the k-th, this may 
serve as a recursion relation. Let us denote by f <K) any solution for 
in this equation — a solution which, in general, will depend on lower 
orders of 5 and on their derivatives. Then one can always add to this solu- 
tion an arbitrary conjugate gradient v X <K) , since the curl operation 
anyway wipes it out. Thus the general solution will have the form 

^(k) = + l (k) (£) (5-2) 

(k) 

with f any particular solution. For instance, one possible form for 
the next two orders turns out to be 


- 

$x l2) ♦ 

* s (1 -v£ (1) 

(5-3) 

S (5) - 

?X (5) 4 


(5-4) 


The vectors f are not unique — given any one set of such vectors, 
we can always form another one by adding to it expressions having the 
form of conjugate gradients. However, if we have any one such set avail- 
able, all near-identity transformations can be characterised by it. 

For lower orders, one can construct solutions by vector manipulation, 
but this quickly becomes difficult. However, since any direct canonical 
transformation must have the form (5-2), any method for obtaining such a 
transformation must lead to a choice of f . 



- 7 - 




Fortunately, there exist at least two such methods. 


First of all, there's the Lie transform. If one expands it and derives 
a formula for , one finds 


.(k) 

S 



(expression involving lower orders) 


(k) 

Obviously, - Vr ; is our 
for f u ’ o 



and the remaining terms give one choice 


Alternatively , one can untangle the transformation based on (T (P, q) and 
express it as a direct transformation. The result is 


= - V Cf^ f (expression involving lower orders) 

(k) 

Again, this gives a choice for f' and, as it turns out, this form is 
simpler to derive and shorter than the other one. Any of these forms can 
be substituted in the Poincare-Von Zeipel method, giving 


H(y) = H* ( y + 21 ^ ] ) 

(k) 

Depending on our choice of f , we get from this the various orders 
either of the "conventional" or the "Lie - type" generating function 
that transforms our variables to a set in which the Hamiltonian does not 
contain the angle variable. 

(incidentally, if we want to use the inverse expansion 2-2, we must 
substitute in H , not in H ) 

The method can also be extended to situations with a slow explicit 
dependence on time ; the constant transformed action variable then goes by 
the name of adiabatic invariant . Since in such cases we need the "conven- 
tional" generating function for deriving the new Hamiltonian, we must use 
the direct transformation based on CT and not, say, that derived from Lie 
transforms. 

In conclusion, let me say that I suspect that this approach to near- 
identity canonical transformations offers more convenience than others. 

For those interested, I have copies of an article on the subject that just 
appeared in J. Kith, Physics. 

What else it means, I don't know, though it seems that the concept of 
conjugate coordinates has something behind it — perhaps some canonical 
formalism involving complex numbers. The use of the conjugate gradient and 
curl also suggest an interesting geometry in phase space, and J hope all 
this will be investigated in due course. 


(end) 
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Afterthought.* on the Application of Direct Canonical Transformations 


(not presented at the meeting) 




The new forms of canonical transformations — those described here and 
Lie's — provide alternative possibilities for new perturbation schemes 
replacing the traditional method of Poincare and Von Zeipel, While all these 
methods lead to equivalent results, I believe — for reasons outlined below — 
that the traditional method still offers the most convenient approach* 


Any of these methods derives a near-identity canonical transformation, 
characterized in each order by the appropriate order of some "generating 
function" ( <r<"> of (2-5), W (K> of (2-5), of (5-2), etc.) y and 

the calculation usually centers around the derivation, order by order, of this 
function and of the transformed Hamiltonian H* . If the transformation is a 
direct one, the task ends here; in the traditional method, however, an inver- 
sion must still be performed to untangle the relationship between old and new 
variables and to bring it to a "direct" ferm. 


Suppose we used a "direct" transformation 


2 * £ eV k) (*> 


with 


5 


(k) 


vX (k) + I 


(y) 

(A-l) 

(k) 

(A-2) 


Then, as on page 7, the basic relation is 


H(y) = 


H*(z) 

“*(? + [?X (k) + £ (k) ]) (*-5) 


and this must now be expanded order by order as in the ordinary Poincaref- 
Von Zeipel method. Different choices of f m <K) lead to different perturbation 
methods: for instance, choosing as in (5-5) 


f^ « 7 ^. V £ (l) 


U-4) 


( vanishes) gives something equivalent to the Lie transform method. On 

the other hand, let the notation 



«. n _ 


V (k) 

S 3 ■ 

denote the particular direct 
transformation generated by >T( P, q) of (l-3), after it is untangled (this 
is the transformation mentioned as the 2nd alternative on top of page ?)• With 
this choice 



Comparison shows that (A-6) has only half as many differentiations as (A-4) 
and therefore ought to be easier to use. Similar simplifications exist in 
nigher orders and for that reason, of all the direct transformations of form 
(A-2), that one which is given in (A-5) appears to be the most convenient one 
for use in a perturbation theory. 

If instead of all this we used the traditional Poincare-Von Zeipei method, 
we would have expanded (l-6) and derived from it CT VK, (£» <j) and H* order 
by order. Having derived q)t we would than invert the relation and 

derive the corresponding direct transformation, simply by writing p in place 
of P everywhere in the generating function and substituting the functions 
(TXp, q) thus obtained in (A-5). 

This procedure derives the same function (within the freedom allowed 
by the method) as does (a- 3) with the choice (A-5) for the direct trans- 
formation — which, as has been pointed out, is probably the most convenient 
choice for (A-3). The difference is that in the traditional method, G -1 *) is 
first derived to all ordeis and then ? tM is formed and the direct transfor- 
mation (A-5) is obtained. By contrast, 'in the direct form of the method (A-3), 
the vectors ■? t *’ (taking the role of ) are already introduced in the 

expansion of the Hamiltonian and have to be carried along throughout the 
entire calculation. 

One may thus conclude that the traditional method is more convenient . since 
it introduces only in the end. Of course, the difference only appears for 

k 2, because f vanishes. 


VCT (k \p, q) + (A-5) 

transformation corresponding to the "conventional" 



The formulas for deriving 


are listed the following page. 
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Note on '^he use of ( A— 5) 


(k) 

The general formula for i ' la (see "Direct Canonical Transformations" , 
to appear in J. Kath. Phys. ) 

k-1 

(k) = - £, s (l) . v cr (k -“ ) (A-7) 

m=l 

(m) 

where S’ - are expansion operators (the asterisk denotes their 
operation), with the fir3t lew being 



o(l) 

= 

V 




s (2) 

=5 

Jl^V + 

i JT^Jl VV 

(A-8) 


s (3) 

= 

71 ^V' + 

71^71^* VV' 





+- 

(1/6) j^( 1 )jj( 1 )^( 1 ) j . > 

V V v> 

and so 
vector 

forth, and 
$'-*>, that 

idle re 
is 

( 8 ) 

71 is just the momentum-like part of the 


J1 ^ 

= ( 

> (p) 

bl * * •• 

il’K o, ... 0) 

(A-9) 


For instance, for the order following (A-6) one gets ( 

^(3) _ xt ( 2 ^ ^d) * 7 

- (J l (l )v)( (A-10) 

This can be expressed directly in terms or the 0"^* : alone, but 

the result ic somewhat more cumbersome than the preceding formula. 



